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RESEARCH MEMORANDUM 


FLUTTER CHARACTERISTICS OF SWEPT WINGS 
AT TRANSONIC SPEEDS 
By Laurence K. Loftin, Jr. 


SUMMARY 

An experimental study of the flutter characteristics of swept 
wings is being mad.e in the Langley transonic blowdown tunnel. The 
purpose of this work is to determine the important effects of some of 
the plan- form variables and to provide the basis for a means of flutter 
prediction at transonic speeds. The investigations have consisted of 
studies of the effects of variations in sweepback angle, aspect ratio, 
and taper ratio through a Mach number range extending from about 0.8 
to 1.35* The investigations have shown that although some further work 
is no doubt required, the basic effects of these plan-form variables 
are fairly well defined. Variations in the center-of-gravity position 
have been shown to have an important effect on flutter at transonic 
speeds. A method of analysis has been developed which accounts for the 
effect of center-of-gravLty position and which indicates the important 
aerodynamic parameters influencing flutter of a certain class of wings 
at transonic speeds. 


INTRODUCTION 


An experimental study of the flutter characteristics of swept 
wings is being made in the Langley transonic blowdown tunnel. The 
purpose of this work is to determine the important effects of some of 
the plan-form variables and to provide the basis for a means of flutter 
prediction at transonic speeds. The investigations have consisted of 
studies of the effects of variations in sweepback angle, aspect ratio, 
and taper ratio through a Mach number range extending from about 0.8 
to 1.55- The variations in the individiial plan-form parameters were 
chosen with a view toward bracketing the range of practical interest 
and consisted of sweepback angles from 0° to 60 °, aspect ratios from 
2.0 to 6.0, and taper ratios from 1.0 to 0.2. The effects of varia- 
tions in certain wing mass parameters have also been briefly stvidied. 
Some of the results of these investigations have been reported (refs. 1, 
2, and 5 ) f whereas other more recent data are not yet generally 
available. The present paper will attempt to summarize and correlate 
some of the trends shown by the resiilts. The experimental techniques 
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employed in the investigations are fully described in references 2 and 
5 and will not he discussed here. 


SYMBOIS 


A 

A 

A 

M 

'^EXP 

Vref 

Vexp 

Vref' 

^2 


Subscripts: 

CG 

ac 

M 

.8 

N 


sweep angle of wing quarter-chord line 
wing aerodynamic aspect ratio 
wing taper ratio 
lifech number 

ratio of experimental to calculated flutter speeds 

ratio of experimental to modified calculated flutter speeds 


ratio of measured coupled second bending frequency tq 
first torsion frequency 


lift-curve slope 


distance along wing chord measured from leading edge, 
fraction of chord 

wing chord length normal to quarter-chord line 


center-of-gravity position 
aerodynamic-center position 
Mach number 

stream Ifech number of 0.8 

direction normal to quarter-chord line 
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METHOD OF ANALYSIS 


Before discussing the trends shown hy the results for the various 
wings, a few remarks regarding the method of presentation and the 
definition of flutter speed may be appropriate. Some typical results 
for a swept wing aare shown in figure 1. A definition of the exact plan 
form is not important in this case. The ffech number is plotted along 
the abscissa and the ratio of the experimental to a calculated, or 
reference, flutter speed is on the ordinate. The reference flutter 
speeds were determined from a Rayleigh type of analysis in which the 
flutter mode was represented by the superposition of the uncoupled 
modes of a cantilever beam and in which the aerodynamic coefficients 
were two-dimensional, incompressible values taken normal to the quarter- 
chord line (ref. 4 ). The necessity of employing such a normalizing 
factor as Vrep in the presentation of experimental flutter results 
seems unavoidable because of the large number of mass, elastic, geo- 
metric, and aerodynamic variables involved. Thus, by. use of a reference 
flutter speed, the mass and stiffness properties of the models and the 
air density, all of which have a profound effect on the actual flutter 
speed, do not appear explicitly in the comparison of the various wings 

Vpyp 

but are implicit in the values of Vp^Ep. Curves of against 

VrEF 

Iifech mmiber, therefore, show the departure of the actual flutter speed 
from a known reference level as a function of Mach number. 

In many of our tests, difficulty is experienced in selecting a 
unique boundary which separates a condition of flutter from a condition 
of no flutter. The data points through which the solid line is faired 
in figure 1 indicate a condition of continuous flutter. The cross- 
hatched area represents a region of doubt in which the behavior of the 
model is characterized by random oscillations and intermittent bursts 
of flutter. As can be seen, the region of intermittent flutter is 
primarily associated with the supersonic range, although this is not 
always the case. The significance of this region of doubt and the 
extent to which the bursts of intermittent flutter may be due to 
excitation by tunnel turbulence in a region of low, but not zero, aero- 
dynamic damping are open to some question. In any case, the flutter 
boundaries to be presented in succeeding figures correspond to the 
condition of continuous flutter as Illustrated by the solid line in 
figure 1. 


RESUITS AND DISCUSSION 


Tbe results of the study of the effect of varying sweep angle are 


shown in figure 2 where the flutter- speed ratio 


^EXP 

Vref 


is plotted against 
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Mach nmber. The wings had an aspect ratio of 4.0, a taper ratio of 
0.6, and were about 4 percent thick in the streamwise direction. The 
sweep angles are seen to be 0°, 50°, 45°, 52.5°^ and 6o°. In the Mach 
nvunber range below about 0.9^ the agreement between experimental and 
calculated flutter speeds is very good, in spite of the oversin5)lified 
representation of the aerodynamic forces in the calculations. As indi- 
cated in the key of the figure, only two modes were employed in the 
calculations for the wings of 0° and 50° sweep. These were the 
\mcoupled first torsion and first bending models of a vtniform cantilever 
beam. In addition to these modes, the second uncoupled bending mode 
\ras employed in the calculations for the other wings. The necessity 
for employing a third mode was found to be closely connected with the 
value of the ratio of second bending to first torsion freq\aency. The 
frequencies forming this ratio were the measured coupled values. The 
third mode appeared to be necessary in order to obtain good agreement 
between calculated and experimental flutter speeds when the ratio of 
second bending to first torsion frequency was in the vicinity of, or 
below, 1.0. 

For Mach numbers greater than about 0.9> the value of the flutter- 
speed ratio increases with Mach number by an amount which depends on 
the sweep ajigle. Very little increase is noted for the 6o° swept wing, 
with progressively more increase accompanying decreases in the sweep 
angle from 6o° to 50°. An inversion in this trend is noted in the 
curve for the unswept wing which falls below the curves for the 50° and 
45° swept wings. No entirely convincing reason for this behavior is 
apparent at the present time, although one possibility suggests itself. 
Difficulties with static divergence were experienced with some of the 
linswept-wing modes. These divergent tendencies may have, in some way, 
obscured the true zero-angle-of -attack flutter boiandary. In any case, 
investigations are now being made of wings of about 10° sweepback angle 
in an effort to clarify these results. 

Some effects of aspect ratio are shown in fig\ire 5 in which the 
^EXP 

flutter- speed ratio is again plotted as a function of Mach number. 

The data shown are for 45° sweptback wings having aspect ratios of 2, 

4, and 6. The taper ratio is 0.6 and the airfoils are 4 percent thick. 
The calciilated and experimental flutter speeds agree quite well at 
subsonic speeds for the wings having aspect ratios of 4 and 6. The 
calculated flutter speeds for the wing of aspect ratio 2 are> however, 
considerably lower than the experimental values at subsonic Mach numbers. 
The discrepancy between experimental and calculated flutter speeds in 
this case is perhaps due to the inadequacy of the two-dimensional aero- 
dynamic coefficients employed in the calculation. At supersonic speeds, 
all three wings are characterized by values of the flutter-speed ratio 
which increased with increasing Ifech number. The shape of the ciorves 
is, however, somewhat different for the three wings. 
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Some indication of the effect of taper ratio is provided in 
figure 4 in which the flutter-speed ratio is plotted against Mach 
number for 45° sveptback wings having taper ratios of 1.0, 0.6, and 
0.2. The aspect ratio was 4.0 for all three wings and the airfoils 
were 4 percent thick. At subsonic tfech nvimbers, the agreement bet\reen 
calculated and experimental flutter speeds is seen to be good for the 
wings with taper ratios of 0.6 and 1.0. The calculations for the wing 
with a taper ratio of 0.2, however, give a flutter speed which is too 
low by about 20 percent. T)ie flutter mode for these wings was charac- 
terized by high frequencies, between the still-air coupled second 
bending and torsion values, with large tip deflections. The first and 
second xmcoupled bending and first vmcoupled torsion mode shapes of a 
beam with a taper ratio of 0.2 were employed in the calculations. The 
flutter mode shape, however, may not have been adequately represented. 
Also, the still-air vibration modes for the wings having a taper ratio 
of 0.2 were highly coupled, which raises some question as to the approx- 
imate method employed for deducing the uncoupled torsion frequency from 
the coupled values. Consequently, the subsonic level of the 
Vryp 

curve for the plan form with a taper ratio of 0.2 is not too well 


VreF 

established 


The data show, however, that the 


Vexp 


curve tends to 


Vref 

rise more steeply with Increasing Mach number as the taper ratio 
decreases . 


The results presented in the preceding three figures all show an 

Vryp 

increase in the flutter- speed ratio — as the Mach number increases 

''REF 

into the supersonic range. The fact that the agreement between calcu- 
lated and experimental flutter speeds becomes poorer as the Mach nxanber 
increases is not svirprising because no account was taken in the Vref 
calculations of the effects of compressibility on the aerodynamic char- 
acteristics. The changes in aerodynamic characteristics with tfech 
number would seem to be primarily a function of wing-plan-form shape. 

An ingjortant question arises, however, as to whether the curves of 

Vexp 

y " . against Mach number are a fvmction only of wing plan form or 

whether these curves may be altered by variations in some of the mass 
and elastic properties of the wing which are hidden in the Vrrp calcu- 
lation. Some understanding of the lii5)ortant aerodynamic parameters 
affecting the flutter speed may be obtained from the simple flutter 
formula given by Theodorsen and Garrick in reference 5* This empirical 
formula is based on the results of low-speed studies of two-dimensional 
wing flutter and is applicable to cases in which the ratio of first 
bending to first torsion frequency is small. A consideration of only 
those elements of the formula which contain the aerodynamic character- 
istics of the wing indicates the following iuiportant proportionality: 
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where the symbols have the following meaning: 


1 


(f) 

CO - (f) 

ac 


V 

% 

(t) 

(f) 


CG 


flutter speed 
lift-ciarve slope 

section center-of-gravity position 
section aerodynamic -center position 


( 1 ) 


The assumption is now made that the departure of the curves of 


'EXP 


from 1.0 as the tfeich nvmiber increases is a function only of the well- 
known rearward shift in the aerodynamic center and reduction in lift- 
curve slope. On the basis of this ass\JH 5 )tion and with the use of the 
relation (l), the following expression for the flutter- speed ratio is 
obtained: 


Vref 


vexp I 

^Iu).8 [©CO - 

1 -8 

Vref 

f(fV -(2) 


11 

^\C/cG \C/a( 



( 2 ) 


The subscript .8 refers to a Nfe.ch nimiber of 0.8 for which the flutter- 

speed ratio . J! 2 ^ . is visually about 1.0 and the corresponding aerodynamic 
Vref 

center is near the 25 -percent- chord station assumed in the Vp^F calcu- 
lation. The subscript M refers to some Mach mmiber higher than 0.8. 
Equation ( 2 ) clearly shows that the reduction in lift-curve slope and 
rearward movement of the aerodynamic center which accompany an Increase 
In Mach number beyond 1.0 should cause an Increase In the flutter- speed 
Vgxp 

Equation (2) also shows that the magnitude of the effect 


ratio 


Vref‘ 


of rearward movements in the aerodynamic' center on the flutter- speed 
ratio depends upon the position of the section center of gravity. 


In order to obtain some indication of the correctness of these 
ideas, a short experimental investigation was made of three wings having 
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identical plan forms but different section center-of -gravity positions. 
The wings had a sweep angle of 45°, an aspect ratio of 4.0, a taper 
ratio of 0.6, suid 4-percent-thick airfoil sections. The wings had 
section center-of -gravity positions of 34, 45, and 57 percent chord. 
TTip results of the investigation are shown in figure 5 in which the 

^EXP 

flutter-speed ratio is plotted as a function of Mach number. 

VREF 


For any given supersonic Mach number^ the value of 


^EXP 

^REF 


is seen to 


increase with forward movements of the center-of -gravity position. In 
fact, the wing with the most forward center-of-gravity position could 
not be fluttered at all, within the opera,ting limits of the tunnel, 
above a Mach number of approximately 1.2. No-flutter points for this 

vp 

wing are indicated by solid symbols. The higher values of for 


the more forward center-of-gravity locations are entirely consistent 
with equation (2) . 


Equation (2) suggests certain possibilities for generalizing the 
data of figure 5 to Include other center-of-gravity positions. The 
values of the lift-ciu*ve-slope ratio and the aerodynamic-center positions 
appearing in relation (2) are imknown and must be fo\md. One possibility 
is to use overall wing lift-curve slopes and aerodynamic-center positions 
as determined from static aerodynamic tests of rigid wings. Such a pro- 
cedure does not yield good results, however, because the deflection of 
the wing is not considered. Another possibility is to regard again the 
aerodynamic parameters appearing in equation (2) as liaraped or integrated 
values and to determine these values with the use of the flutter data of 
figiare 5 and equation (2). This procedure has been followed herein. 

The aerodynamic-center position and the lift-curve-slope ratio at any 
given Mach number are assumed to be a function only of the plan form 
and, hence, would be the same for the three 45° swept wings having 
different center-of-gravity positions. This assimiption implies that 
the flutter mode shapes for the three wings are not markedly different. 

Equation (2) indicates that the difference in the curves of at 

VREF 

supersonic speeds for the three AvLngs of figure 5 is expressed by the 
difference in value of the ratio 



8 


CGHFiDElSTtAL 


••• ^ •• 
• • • 

• • • • 

• • • 


MCA RM L55E19a 


whereas the lift-curve-slope ratio 


(^^)m 


has the same effect on ^all three wings • On the hasis of this synthesis 
of the effects of lift-curve slope and aerodynamic center, the faired 
ciu:ves of figure 5 for the three wings were cross -plotted in such a way 
as to determine the variation of the lumped, or effective, values of 
the aerodynamic-center position and lift-curve-slope ratio with Mach 
number. The resulting aerodynamic-center position is shown as a 
function of Mach number normal to the quarter-chord line in figiare 6. 

The aerodynamic center is seen to shift from the 0.25c station to the 
0.54c station as the normal ^fa,ch number varies from 0.55 to 0.95- These 
values appear qviite reasonable. The variation of the Uft-curve-slope 


ratio 


(^ig).8 


with Mach number normal to the quarter- chord line is 


shown in fig\u*e 7* The ratio of the lift-curve slopes at stream Mach 
niimbers of 0.9 and 1.2, as determined from some unpublished static 
aerodynamic tests of a rigid 45° sweptback aspect-ratio-4.0 wing, is 
shown by the symbol in this figure. 


An indication of how well the deduced variations of aerodynamic- 
center position and lift-curve ratio describe the results obtained for 
the wings with different center-of -gravity positions is provided in 

figure 8. In this figure, the ratio is plotted against stream 

^REF 

Mach number for the wings with different center-of -gravity positions. 
The values of Vj^' were determined from eqviation (2) by letting 


Yref * _ 1 


[(I)CG - ' 

0.25 

.8 

VreF 

(%)„ 

(!) 

L. ■ 1 

CG 

(-) ' 
'‘C'ac 

. M 


(5) 


and using the values of 



and 




as given in figures 6 


and 7 . The correlation is excellent, with no systemtic trends evident 
for the wings with different center-of -gravity positions. 


The variations of aerodynamic-center position and lift-cxirve -slope 
ratio shown in figures 6 and 7 , respectively, were determined from 
flutter tests of a particular wing plan form. Application of the results 
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to an arbitrary wing plan form is not, in general, permissible. For 
the restricted case in which the sweep angle is the only plan-form 
parameter varied, however, one might expect the values of aerodynamic- 
center position and lift-curve-slope ratio for one sweep angle to be 
roughly applicable to other sweep angles on the basis of equal Mach 
members normal to the quarter-chord line. On the basis of this rather 
crude assumption, the aerodynamic-center position and lift-curve-slope 
data of figures 6 and "J, together with the formula (5)> have been used 
in an attempt to correlate the swept-wing data of figiire 2. The results 

are presented in figure 9 in the form as a function of Mach 

Vref' 

number normal to the quarter-chord line. Also included in figinre 9 are 
the data for the 45° swept wings with center-of -gravity positions of 54 
and 57 percent chord. The unswept-wing data of figure 2 are not 
included. For normal ^fe,ch mambers less than 0.7^ the correlation is 
within the scatter of the data for individual wings. At higher Mach 
numbers, the correlation is not quite so good, with the maximum dis- 
parity between the data points and the line = 1.0 being about 

Vref* 


15 percent. 


The correlation of figure 9 indicates that, at least for the class 
of wings considered, the aerodynamic center and lift-curve slope are 
the important aerodynamic characteristics controlling the variation 

of with Mach number. The fact that the correlation, figure 9^ 

Vref 

was achieved without any consideration of the effect of compressibility 
on aerodynamic lag is perhaps also of some significance. The formula, 
equation (5), together with the values of aerodynamic-center position 
and lift-curve -slope ratio given in figures 6 and J, respectively, may 
perhaps prove of some use in estimating the effect of variations in 
center-of -gravity position on the flutter speed of wing plan forms of 
the same general class as those considered. The generality of the 
method, in an absolute sense, is however difficult to access. 


CONCLUDING REMARKS 


Transonic flutter investigations have been made of swept wings 
having different sweep angles, aspect ratios, and taper ratios. 

Although some further work is no doubt required, the basic effects of 
these plan-form variables seem fairly well defined. Variations in the 
section center-of -gravity position have been shown to have an inportant 
effect on flutter at transonic speeds. A method of analysis has been 
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developed which accounts for the effect of center-of-gravity position 
and which indicates the in^jortant aerodynamic parameters influencing 
flutter of a certain class of wings at transonic speeds. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., April 27, 1955, 
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FLUTTER BOUNDARIES 



0 .8 1.0 1.2 1.4 

M 

Flgiire 1 





NACA RM L55E19a 




• •• •• 


EFFECT OF TAPER 

A. 45 '; A* 4.0 



L I I I I 

0 .8 1.0 .. 1.2 1.4 


Figure if 


Vexp 

Vref 


EFFECT OF CENTER OF GRAVITY 
A»45*, A»4.0j x=0.6 
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Figure 5 
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AERODYNAMIC -CENTER POSITION 
A*45*j A *4.0; X«0.6 



Clq ratio 

A»45'; A = 4.0; X«0.6 



Mn 


Figure 7 
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CENTER-OF-GRAVITY CORRELATION 

A*45®j A=4.0i X=0.6 




Figure 8 


SWEPT-WING CORRELATION 
A. 4.0; X*0.6 
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Figure 9 
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